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Question 1 : Quantification 

Expliquer le principe de la quantification. Définir les termes suivants : pas de quantification, niveau de 
représentation, seuil de décision, saturation d'un quantificateur. Comment calcule-t-on le SNR en sortie 
d'un quantificateur uniforme ? Expliquer la démarche de conception d'un quantificateur optimal. 
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La quantification est le procédé qui permet d'approcher un signal continu par les valeurs d'un 
ensemble discret. La quantification correspondant à des niveaux discrets, il en résulte une perte 
d'information, pour peu que l'on puisse mesurer des variations inférieures aux longueurs 
d'intervalle de quantification avec des moyens analogiques. 

Nous allons considérer que la quantification est instantanée et sans mémoire, ce qui veut dire 
que les échantillons pris en un temps "t" sont indépendant de la valeur des échantillons 
précédents et à venir. 

Les seuils de décision forment l’ensemble de valeurs de représentations possibles. 


Exemple : 

Je numérise avec un pas de 2V un signal dont l’amplitude va de 0V à 10V . Mes seuils de 
décisions seront 0-2-4-6-8-10 V. 


Les niveaux de représentations sont l’ensemble d’amplitude discrète du signal quantifier. 

Le pas de discrétisation est la distance entre deux niveaux de représentation, autrement dit la 
plus petite variations du signal que le système peut distinguer. 

Le seuil de décision du quantificateur est la valeur limite pour laquelle le signal au temps t est 
quantifier à une valeur supérieure ou inférieure 

La saturation d’un quantificateur survient lorsque le signal d’entrée excède la plage d’entrée du 

quantificateur. Une fois que la plage du quantificateur a été dépassée, la sortie du quantificateur 
restera à la valeur minimale ou maximale jusqu’à ce que l’entrée revienne dans la plage 































SNR Signal Noise Ratio 

Soit une variable aléatoire continue M compris entre [~m max , m max ], avec un quantificateur uniforme, 
le pas de quantification est 

Où L est le nombre de niveau de représentation du quantificateur 

L'erreur Q considéré comme Q = M - V est bornée entre < Q < § et si A est assez petite, 
l’erreur Q peut être considérée comme une variable uniformément distribuée. La densité de 
probabilité de Q est alors 


f Q (q) = A P° Ur ~ 2 ^ 1 ^ 2 

Ceci est vrai si le signal ne sature pas le quantificateur 
Avec une moyenne de quantification nulle, la variance vaut 
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Si L est le nombre de niveau de représentation, le nombre de bit par échantillon sera de 2 R (R le 
nombre de bits utilisés par échantillon) et donc on peut écrire 


A = 



La variance peut donc être réécrite selon cette forme 
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Alors si P est la puissance du signal correspondant au message, on peut exprimer le SNR du 
quantificateur comme 


SNR 0 
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On observe donc que le SNR 0 du signal de sortie du quantificateur croît exponentiellement avec le 
avec le nombre de bit par échantillon R. 





Design d’un quantificateur 

"Sélectionner les niveaux de représentation et les intervalles de partition de manière à minimiser le bruit 
de quantification pour un nombre de L donné." 

On cherche donc un ensemble de "v k " et de "J k " qui minimisent la distorsion moyenne. 

D = ^2 d(m,v k )f M {m)dm 

k =1 ' meJ k 

où d(m,Vk) est la fonction de distorsion = (m - v k ) 2 

Un quantificateur est en réalité composé de 2 partie, un encodeur (émetteur) et un décodeur 
(récepteur). 

Option 7 ; Design de l'encodeur pour un émetteur donné 

On suppose donné un "Codebook” regroupant l’ensemble des "v k " utilisé pour le décodage voulant 
minimiser la distorsion moyenne, 

D= I d{m,g(m))f M {m)dm>y' [ [min d(m, v k )]f M {m)dm 

J-A ^ J me J k v k eC 


Option 2 : Design du décodeur pour un encodeur donné 

Connaissant les intervalles, on peut trouver un codebook optimal minimisant la distorsion moyenne: 

D = J2 (m - v k ) 2 f M (m)dm 

/c=1 J 


L'optimum se trouve là où la dérivée de cette fonction vaut 0: 

= -2 V / (m - v k )f M {m)dm = 0 

° Vk k =1 Jmej k 



I '» -f M t m ) dm 

mcJ t _ 




où le dénominateur est la probabilité que l'échantillon appartiennent à l’interval J Â , 







• Additif car il est ajouté au bruit intrinsèque du système. 

• Blanc car sa puissance est uniforme sur toute la bande fréquentielle du 
système. C'est une analogie avec la couleur blanche, qui correspond une 
émission uniforme sur toutes les fréquences du spectre visible. 

• Gaussien car il a une distribution normale dans le domaine temporel avec 
une moyenne nulle (selon la courbe de Gauss). 

Question 2 : Filtre adapté 

Expliquer l'intérêt de la notion de filtre adapté dans la conception d'un système de communications 
numériques. Calculer la probabilité d’erreur binaire obtenue en sortie du filtre adapté, pour un canal 
AWGN. Additive White Gaussien Noise, bruit additif blanc Gaussien 
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Objectif du filtre adapté : récupérer un signal corrompu par du bruit blanc (Densité Spectrale de 
puissance = constante) fréquent en Télécom 

Soit le signal g(t) additionner de bruit blanc w(t) correspondant au signal résultant x(t) en sortie du 
canal 
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x(t ) = g{t) + w(t ) avec 0 < t < T 

Le filtre h(t) a alors comme but de détecter le signal en t = T suivant l’équation 

y(t) = x(t) ® h(t) = g(t) + n(t) 

Attention, le bruit résultant du filtrage n’est plus blanc 

On peut alors définir le SNR à maximiser pour que le bruit vienne perturber le moins possible le signal 
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Avec g 0 (t ) le signal reçu (déterministe) et le bruit n(t) qui est aléatoire 


Numérateur de l’expression du SNR 
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Dénominateur de l’expression du SNR 
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Inégalité de Schwarz 

L’inégalité suivante est toujours vérifiée pour deux fonctions complexes d’énergie finie 
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Avec l’égalité dans les cas où k est constant (arbitraire) 

Cpj (x) = k(p* 2 (x) 











Nous retrouvons donc l’expression du SNR 
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Le SNR est donc indépendant de H(f) (le filtre). On peut ensuite trouver la réponse optimale du filtre 
pour que l'inégalité de Schwarz devienne égale 
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Propriété : le ffiltiii signal ft'inrit n en sortie dépend que du ffppgtf entre 

l’énergie E du signal d’entrée et la densité spectrale de puissance N 0 du bruit blanc en entrée 

Spectre du signal en sortie 

G 0 (f) = H opt (f) .G(f) = k \G(f)\ 2 e-^ fT 
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Energie du signal 
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Puissance moyenne du bruit en sortie 
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SNR final 
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Probabilité d'erreur due au bruit 

Prenons un train d'impulsions binaire x(t) de type "polaire NRZ" où chaque symbole 0 ou 
durée T h et où (+V etV sont les amplitude A respectivement). 


Tension 
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Le signal reçu et affecté par le bruit blanc vaut alors sur un intervalle 0 < T b < t : 


x(t) = A + w(t) symbole(l) 
x(t) = A + w(t) symbole(O) 

Supposons que le récepteur est parfaitement synchronisé, c'est-à-dire qu'il connaît parfaitement la 
position des symboles ainsi que leur durée T b . Si ne n'est pas le cas, il y aurait erreur sur 
l'interprétation à la sortie. En effet, le canal où l'on envoie le train d’impulsion génère un certain 
déphasage sur le signal envoyé. 

Sur une longue séquence de symboles identiques, le récepteur doit être synchronisé avec 
l'émetteur afin de compter le bon nombre de symboles ! Par exemple, si un émetteur envoie 10 bits à 
0 consécutivement et que la cadence de l'émetteur est de 400Hz (soit T h - 2.5ms), la séquence de 
10 symboles durera 25ms. Si le récepteur est parfaitement synchronisé, il comptera un symbole 
toutes les 2.5ms et atteindra également 10 symboles échantillonnés. Si par contre le récepteur n'est 
pas synchronisé et considère par exemple la durée d'intervalle T = 2ms il comptabilisera 
respectivement 12 symboles au lieu de 10, ce qui fausse le résultat ! 

Supposons également que la forme des impulsions est connue (impulsions rectangulaires 
d'amplitude A) par le récepteur. L'objectif est de savoir déterminer si le symbole est 0 ou 1 sur base 
du signal bruité observé par le récepteur. D'où l'apparition de notion de probabilité d'erreur due au 
bruit et du taux d'erreur binaire. 


PCM wave 
s(») 



Le signal s(t) de type PCM (Puise Code Modulation) est l'entrée de notre système. Ce type de 
modulation est utilisé pour échantillonner la voix en télécommunications et pour le son. Sur base du 
schéma, on voit que le signal PCM (s(t)) est bruité par du bruit blanc Gaussien w(t). Au travers d'un 
filtre adapté, on échantillonne le signal réceptionné (aux temps T b ) qu'on analyse en le comparant à 
une valeur de seuil K. Si la valeur est supérieure à À, alors on considère que le symbole échantillonné 
vaut 1, sinon vaut 0. 




























Ce seuil À n’est pas forcément juste, et on distingue deux types d'erreur lorsque le symbole de sortie 
est erroné : 

Erreur de type 1 : choix du symbole 1 au lieu de 0 
Erreur de type 2 : choix du symbole 0 au lieu de 1. 

Pour la transmission d'un symbole 0, à la sortie du filtre adapté on obtient le signal y(t) échantillonné 
sur base x(t) : 

T .b tu 

y (0 = J x (0 dt = -A + J w(t) dt 

Th o r ‘o 

Le deuxième terme de l'addition est en moyenne nulle car pour rappel l'AWGN induit une densité 
spectrale de puissance constante. 

Le signal y(t) faisant office de variable aléatoire Gaussienne de moyenne -A et de variance (espérance 
mathématique de la puissance du signal divisée par la longueur d'intervalle au carré) : 

V Ÿ V Ÿ V ? 

o v 2 = y ~ 2 E [ J J xv (t) w(u) dt du\ = J J E [ w ( t ) w{u)\ dt du = ^ J J R w (t, u) dt du 
b o o 6 o o b o o 

L'espérance des sommes étant la somme des espérances et avec E [w(t) w(u)] = R w (t, u) étant la 
fonction de covariance à l'instant t et à l'instant u. Le bruit blanc ayant une PSD = N 0 , sa fonction de 
covariance vaut : 

R w (t, u) = N -fUt-u) 

La covariance du bruit blanc est donc une impulsion de Dirac et la variance de y devient donc : 



Voyons maintenant la fonction de distribution conditionnelle de y : 
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En analysant les gaussienne, on distingue : 

p l0 , la probabilité de décider 1 alors qu'on a transmis un 0 
p 0l la probabilité de décider 0 alors qu'on a transmis un 1 





On espère que ces probabilité soient la plus petite possible. Il faudrait donc que dans le cas où on 
envoie beaucoup plus de 0 que de 1 on ai A le plus grand possible pour minimiser p l0 et K le petit 
possible pour minimiser p 0l . 

Le taux d'erreur binaire est le rapport du nombre d'erreurs binaires (nombre de symboles erronés) au 
nombre total de symboles transmis. Par exemple, si les 10 symboles 1001101110 sont envoyés, et 
qu'en sortie le récepteur détermine la trame 1011001100, on voit que 3 erreurs 
ont été commises et le taux d'erreur binaire vaut alors 30%. 

Le pire taux d'erreur binaire est 0.5 (50%). En effet, si tous les bits sont erronés, alors létaux d'erreur 
binaire vaut 1 (100%) mais il suffit d'inverser tous les bits et n’y a plus la moindre erreur. 

En mathématique, la fonction d'erreur (aussi appelée fonction d'erreur de Gauss) est une fonction 
entière utilisée en analyse. Son expression est la suivante : 

erf(x) =zje~ ,2 dt 
o 

Cette fonction erf(x) est symétrique par rapport à l'origine des axes et son domaine est compris entre 
-1 et +1. La fonction complémentaire de cette fonction d'erreur est : 
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Probabilité d'erreur de type 1 

A savoir que y >Â alors qu'un 0 a été transmis. Avec l’expression de la gaussienne f (y | 0 ) on 
calcule : 


GO 

P 10 = \fy(y I 0 )dy 


■\jnN 0 /T b 


dy 


î 

r o/T„ 



dy 


Pour faciliter l'expression on change de variable en posant 
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Probabilité d'erreur de type 2 

A savoir que y < À alors qu’un 1 a été transmis. On obtient de manière similaire : 


a ™ = °- 5 er f c (-wrr) 

01 V-^o 1 1 b 

Ainsi, on peut calculer la probabilité d'erreur moyenne Pe : 

p c = PoP 10 + PlPoi 

Il s’agit de la probabilité d'être en présence d'un bit 0 et de choisir 1 additionnée à la 
probabilité d'être en présence d'un bit 1 et de choisir 0. 

P, - 0.5p 0 erfc (-j===) + 0.5 p, 


Le seuil optimal X OPT peut se calculer comme suit : 

Vr = âatI 1 o S(pJP\) 


En effet, on remarque que le cas particulier où p 0 = p l = 0.5, c'est-à-dire dans le cas 
d'équiprobabilité où on a 50% de chance de tomber sur un 0 et 50% de chance de tomber sur un 1, le 
seuil optimal vaut X OPT = 0 car log(1) = 0 


Dans le cas où p m = p l0 on est face à un canal "binaire symétrique". La probabilité moyenne en 
considérant ce cas là et le cas précédent ( p 0 = p x et donc À = 0) vaut : 
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Soit E b l'énergie du signal transmis par bit : 

E b = A 2 T b 


On obtient alors la forme finale de Pe : 


P e = 0.5 erfc (^E h /N 0 ) 

La probabilité d'erreur symbole dans le cas d'un canal binaire symétrique dépend donc 
uniquement du rapport E b /N 0 entre l'énergie du signal transmis par bit (E b ) et la densité 
spectrale de puissance du bruit (N 0 ) 


Comme le montre le graphique semi-logarithmique ci-dessus, plus le rapport augmente, plus 

la probabilité de faire une erreur diminue ! On note ainsi une 
amélioration exponentielle du taux d'erreur en fonction de E b /N o- 
Alors, pour obtenir un taux d'erreur satisfaisant (inférieur à 10 3 ), une 
augmentation modeste du rapport 
chercher un rapport de 20dB par exemple. 


E b /N 0 suffit, pas la peine d'aller 


Avec P e borné tel que P 
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Question 3 : Critère de Nyquist 

Expliquer la notion d'interférence entre symboles (ISI) et son impact sur les performances des systèmes 
de communications numériques. Dériver le critère de Nyquist et donner un exemple de mise en oeuvre. 
Comment faire pour lutter à la fois contre le bruit (filtre adapté) et contre l'interférence entre symboles 
(filtre de Nyquist) sur un canal AWGN? 
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Le principe de la modulation d'impulsions en amplitude réside dans l'encodage du signal analogique 

par rapport à son amplitude. Le signal est échantillonné à intervalles réguliers et le signal PAM Puise Amplitude 

prendra les valeurs de l'échantillonnage. 

Exemple avec une PAM-4 : un signal sinusoïdal de 3 volts d'amplitude et 4 états (+3V,+1 V,-1 V et -3V) 
donc utilisation de mots de 2 bits pour chaque symbole. 
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Interférence entre symboles 

Soit un système de transmission par modulation PAM de type binaire, c'est-à-dire PAM-2 (un seul bit 
et 2 symboles). Les amplitudes a k des impulsions sont de +1 pour le symbole binaire b k = 1 et -1 
pour le symbole binaire b k = 0. 

La séquence a k passe dans un filtre de mise en forme g(t) pour donner le signal s(t) 

s(t) = I ct k g(t - kT b ) 

k 

Le signal s(t) ne retombe pas à OV à chaque période T h contrairement aux impulsions a k . On arrive 
alors au bout de l'émetteur du système de transmission dont voici le schéma complet : 
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Le signal est donc envoyé vers le récepteur à travers un canal dispersif dans lequel le bruit blanc 
perturbe le signal. Le canal est caractérisé par une réponse h(t) et le bruit blanc par w(t). 

Voyons un peu ce qu'engendre le canal dispersif. Dans le cas précédent on avait considéré un canal 
idéal. La réponse h(t) est maintenant modélisé par un filtre linéaire au canal. Pour le canal idéal la 
réponse est 


h ideaM) = a. 6 (t ~ X) 

Avec a l'atténuation et i le délai de propagation 
En fréquentiel, ce filtre à pour réponse 


H idela (J) 


a. e~ J2nfz 


Le filtre est considéré idéal (pas de déformation du signal) si ces deux conditions sont remplies 


\m\ = « 

d ( arg(H(W) ) 
2 TT ' df 


( distorsion d'amplitude si non respecté ) 
i ( distorsion de phase si non respecté ) 


Au niveau du récepteur, le signal x(t) reçu passe dans un filtre de réception c(t ). Le signal de sortie de 
ce filtre y(t) est alors échantillonné à la même cadence que celle de l'émetteur (f, = iT b ) avec i un 
nombre réel entier (positif). Après chaque période de durée T h on échantillonne la valeur de y(t) et on 
obtient le signal j/(r.). Enfin, les échantillons de j(r ; ) sont soumis à un organe de décision qui les 
compare au seuil A. Exprimons la sortie du filtre de réception y{t) sur base des filtres : 


y(t) = x(t) ® c(t) = x 0 (r) ® c(t) + w(t ) ® c(t) = (i p(t) a k + n{t) 

Avec p(t) résultant de la mise en cascade de 3 filtres (émission g(t) , canal h(t) et récepteur c(t )) 
appliqué au signal d'entrée a k et «(r)un bruit Gaussien filtré (donc plus blanc) 

y(t) = P Z a k p (t - k T b - t 0 ) + n(t ) 

k 

où t 0 représente un temps de propagation (délai) et peut être négligé sans perte de généralité (dû au 
canal). On a aussi p qui est le facteur global de changement d’amplitude, càd que l'amplitude finale 
du signal est p fois plus grande que l'amplitude du signal a k de départ. C'est l'atténuation du signal 
dû au canal. Afin de déterminer p, on normalise l'impulsion p(t) au temps 0 : p( 0) = 1 


Une fois échantillonné, le signal peut s'écrire 


V dû = P% = ! - + I a k p [(/ - k) T b \ + n{t t ) 

k != i 

C'est dans le ième intervalle qu’a été envoyé le bit qui nous intéresse a k=i . Donc ici seul le premier 
terme nous intéresse, les deux autres étant des termes perturbateurs. En l'absence de bruit n(ti) et 
d'interférence entre symboles 

P Z a kP [(' ~ k ) T b\ —> y (h) = f 1 »/ car P ( iT b ~ kT b) = 1 si i = k 

k != i 



La conception d'un système de transmission numérique en bande de base a pour objectif de 
minimiser l'effet du bruit et de l'interférence entre symboles (ISI) de manière à maintenir le taux 
d'erreur le plus bas possible. 

Dans les systèmes à SNR élevé (exemple : le réseau téléphonique), les performances sont 
principalement limitées par l'ISI, et la présence du bruit peut être négligée en première 
approximation 

NB : k est le numéro de l'impulsion PAM envoyé tandis que i est le numéro de l'échantillon de y (r ; ). Si 
k = i cela veut dire que l'échantillon ne correspond pas à l'impulsion envoyée, il y a donc interférence 
avec un autres symbole. 

Critère de Nyquist 

Pour rappel, le récepteur a pour tâche d'extraire la séquence y (; iT b ) en sortie du filtre de réception, 
puis de décoder cette séquence afin de récupérer les symboles binaires b k . Le contrôle du 
canal équivalent p (t) a pour but de satisfaire au mieux la propriété qui suit : 

p(iT b - kT b ) = 1 si i = k car /?(0) = 1 
p (; iT b - kT h ) = 0 si i ! = k 

Afin de n'avoir aucune interférence entre symboles (mais cela ne veut pas dire que cette 
propriété est respectée, il faut tenter de la satisfaire). Si cette propriété est respectée, alors la 
réception sera parfaite (en l’absence de bruit) 

y (*,•) = pa, 

Le critère de Nyquist est simplement l'expression de cette condition de réception parfaite (pas d'ISI) 
dans le domaine fréquentiel. Le spectre de l'impulsion global p(nT b ) échantillonné à cadence jr 

peut s'écrire 



Pour une réception sans ISI, le spectre devient 
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Étant donné que la transformée de Fourier d'une impulsion de Dirac vaut 1. Le critère de Nyquist 
s'écrit finalement 
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On en conclut que si le canal global p(t), composé de la mise en cascade du filtre d'émission, du 
canal physique et du filtre de réception, a un filtre P(f) satisfaisant le critère de Nyquist (ci-dessus) 
alors les échantillons du canal enregistrés à une cadence jr sont exemptés d'interférence entre 
symboles. 

Le canal de Nyquist idéal a un spectre P(f) rectangulaire de bande passante W = fp- 
L'impulsion p(t) correspondante dans le domaine temporel est donné par 

Pif = ” 2 Tm t] = sinc(2Wt) 

Exemple 


Binary sequence 10 110 10 



On remarque bien que les impulsions successives (échantillonnées aux instants t = iT b ,caà t = 0; 
t = T b \t = 2 T b ]... n'interfèrent pas les unes sur les autres car on voit que les 0 sont bien distincts 
par rapport aux 1 et tous les 0 se superposent à intervalle régulier T b . A chaque période ^.toutes 
les courbes passent par 0.0 (ou 1.0), il n'y a donc jamais la moindre confusion possible entre 2 
symboles donc pas d'ISI ! 

Il s'agit du canal à la bande passante la plus faible (fréquence de Nyquist W ). Cependant on constate 
2 inconvénients : 

1. Transition abrupte de la réponse du filtre (rectangle), alors que c’est irréalisable en pratique 

2. L'impulsion p(t) décroît en 1 , ce qui est très lent, dégradant ainsi de façon importante et 
les performances en cas de mauvaise synchronisation. 








Une solution pour palier à ces inconvénients est l'extension de la bande passante depuis la valeur 
minimale W = vers une valeur de B r = [W; 2W] Pour un filtre dont le spectre est limité à la 

bande [-2W ; 2W ], on peut se contenter de spécifier le critère de Nyquist dans la bande [-W ;W] comme 
ceci 


Pif) + Pif ~ 2W) + Pif + 2W) = ^ avec - W </< W 

Une solution répandue à ce problème est le filtre en cosinus surélevé (raised cosine spectrum) 
constitué d'une partie plate et d'une partie de roloff de forme sinusoïdale 


2 WPif) 



£ 

IV 


a le facteur de roloff est toujours compris entre 0 et 1. Il représente la bande passante en excès par 
rapport à W 

a= 1 - - 
u 1 w 

La bande passante totale va alors varier entre W et 2W et vaudra alors 

B t = W ( 1 + a) 

Notes d'observation : Le rectangle contient la bande de canal le plus étroit possible ! C'est celui que 
l'on connaît : le filtre idéal, impossible à avoir en pratique. Mais comme on s’autorise à augmenter la 
bande passante, on observe que : 

Entre -2 à 2 on consomme 100% de plus que la bande la plus étroite 
Entre -3/2 et +3/2 on consomme 50% de plus que la bande la plus étroite 

Donc visuellement, par symétrie si on périodise ces courbes on obtient en additionnant les répliques, 
une horizontale (respect de la condition de Nyquist) 

En règle générale, on a 3 parties dans une courbe dont le facteur roloff a est entre 0 et 1 : 

Un plateau constant 

Un morceau de cosinus de part et d'autre 
Un morceau nulle 

Pour un roloff de 0, les parties sinusoïdale disparaissent tandis que pour un roloff de 1, le plateau 
constant disparaît 


Question 4 : Codes à réponse partielle 

Exposer le principe et l'utilité des codes à réponse partielle. Décrire les systèmes duobinaire et 
duobinaire modifié. Généraliser aux codes à réponse partielle de classe quelconque. Comparer au filtre 
de Nyquist. 
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La technique des codes à réponse partielle ('Correlative-Level Coding' ou 'Partial- Response 
Signaling') permet d'atteindre la limite théorique de 2W symboles par seconde dans un canal de 
largeur W Hz, en utilisant des filtres réalisables (pas de transition abrupte du spectre) et tolérants aux 
erreurs detiming (affaiblissement rapide de la réponse impulsionnelle). 

L'idée est d'introduire volontairement de l'ISI (Interférence entre symboles), de manière 
contrôlée, pour pouvoir l'interpréter correctement dans le récepteur. 

Codage duobinaire 

Le codage duobinaire est un code à réponse partielle de type 1. Il utilise des symboles binaires b k = 1 
ou 0. Pour rappel, la modulation PAM produit des amplitudes a k = +1 ou 1 selon b k . Le codeur 
duobinaire quant à lui va produire des amplitudes c k à trois niveaux : -2, 0 et 2. 

c k = a k + a k-\ 

La séquence d'impulsions décorrélées à 2 niveaux a k est donc transformée en séquence 
d'impulsions à 3 niveaux c k , ce qui revient à introduire de l’ISI de manière contrôlée. La séquence 
d'impulsions c k passe ensuite dans un filtre de Nyquist idéal. 


Input 

two-level 

sequence 



Output 

sequence 

Iql 


Filter W|(/) 


La réponse en fréquence du filtre équivalent donne 

= H nyq (M 1 + e~ 2 ^ =H nyq (f)[e^+ e^]e^ T " = 2H nyq (f) cos (nfT b )e^ 

Pour un filtre de Nyquist idéal (réponse en fréquence rectangulaire) on a donc : 

2 H nyq (f) cos 0 KfT b )e ^ avec \f \ < ^ 


















hj{t) 


T b 2 sin\ÿ\ 
nt(T b -t) 


L'impulsion diminue en d-, ce qui est plus rapide que l'impulsion de Nyquist 



Attention, si le premier bit est erroné, toute la séquence le sera ! Sans précodage, il suffit de soustraire 
a k _j àc H pour obtenir a k 

a k = c k - a k -1 

Cela peut amener des problèmes des systèmes à rétroaction de décision et donc un risque de 
propagation des erreurs - Solution : Précodage (addition modulo 2 (non linéaire)) 

<4 = b k ®d k _ j 


. Table 4.1 Illustrating Example 4.3 on duobinary coding 

Binary sequence (fc t ) 0 0 10 110 

Prcccxlcd séquence {</*) 11100100 

Two-level sequence (<ït) +1 +1 +1 —1 —1 +1 —1 —1 

Duobinary coder output (c*) +2 +2 0 —2 0 0 —2 

Binary sequence obtained by 0 0 10 110 

appiying decision rule of Eq. (4.76) 


Prix à payer : moindre résistance au bruit car plus grand nombre de niveaux d'amplitude 






Par rapport au codage PAM où on avait +1 et -1, on a maintenant 3 états (-2 ; 0 ;+2) 
Statistiquement, la moyenne du signal change étant donnée qu'elle est proportionnelle au carré de 
l'amplitude. 

• Pour la PAM 

On a 50% de chance de +1 et de -1, donc on calcule 


0.5 . P + 0.5 . (-1) 2 = 1 


• Pour le duobinaire 

On a 25% de chance de +2, 50% de 0 et 25% de -2 donc on calcule 

0.25,2 2 + 0.5.0 + 0.25 . (-2) 2 = 2 

On se rend compte que pour obtenir la même chose que la PAM dans le cas duobinaire, il faut diviser 
par V2 (la moyenne devient alors = 1). Dans ce cas, les signaux deviennent comparables. L'écart-type 
est maintenant beaucoup plus faible que (écart de V2) et les spectres se chevauchent plus, ce qui 
laissent à penser que le duobinaire est plus sensible au bruit que la PAM 

Il faut précoder la séquence de bits avant de coder en duobinaire. Ainsi, lorsque c k à 3 états arrive 
dans l'organe de décision, il ne reste plus qu’à décider en sortie b k = 0 si c k = + ou - 2 et 
b k = 1 si c k = 0 (le signe n'a plus d'importance, on retombe de 3 états à 2 états) 



Question 5 : Égalisation 

A quoi sert l'égalisation ? Pourquoi le choix d'un égaliseur MMSE est-il intéressant ? Commenter les 
équations permettant de dériver l'expression de l'égaliseur MMSE. Exposer le principe de l'égalisation 
adaptative et de l'égalisation à rétroaction de décisions. 
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Solution pour le bruit de canal —» filtre adapté 

Solution pour l'interférence entre symbole —* critère de Nyquist 

Sauf que ces deux phénomènes se produisent simultanément, et le filtre de réception doit donc être 
conçu en conséquence 

La solution compromis entre bruit et ISI —*• Minimum Mean Square Error Equalizer (MMSE) 

Échantillon recueilli en sortie du filtre de réception 

oo oo 

y(iT b ) = I a k \ c(x) q(iT b - kT h - x) dx + | c(x)w(iT h - i) ch = Ç. + n. 

k —oo —oo 

Idéalement entre on voudrait recevoir y (; iT b ) = a et donc on définit l'erreur comme 

e i = y(’ T b ) - a i = Ç; + n i - a i 

L'erreur quadratique moyenne est quant à elle égale à 

J = 0.5 E [e t Y = 0.5 E [ Ç. 2 ] + 0.5 E [ n t 2 ] + 0.5 E [ a, 2 ] + E [ l ii n i ] + E [ n i a i ] +E [ Ç ( .a,. ] 

En calculant les termes on trouve 


(1.47) £ [ ^. 2 I = I j cfi,) Rq{x |, i 2 ) c(i 2 ) dx t dx 2 

—OO —OO 

< ? ) < ? ) 

(1.48) £[Ç,. 2 ] =^2 j J c(t,). c(t 2 ) . 5(x j — x 2 )dx l dx 2 

—00 —00 

(1.49) E [a, 2 ] =1 

(1.49) = E[n t a t ] = 0 

00 00 

(1.49) E [ ] = X I = j c(x)q(iT b - kT b - x 5i = j c(x)q(~x)dx 

k -00 —00 

Et donc l'erreur quadratique moyenne est égale à 


00 00 


J = 0.5 + 0.5 J J (R q {t~x) + 0.5N o ô (t-x ))c(t)c(x)dtdx 


—OO —00 


J c(t) q(- t) dt 


Pour trouver le filtre de réception optimal (MMSE), on différencie et on égale à 0 J par rapport a c(t) 


j (R q (t-x) + 0.5N o ô (r — t ) ) c(x) dx = q(-t ) 


Solution dans le domaine fréquentiel 

(S q (j) + N 0 /2) C(f) = Q* if) 

Ce qui donne une solution explicite pour le filtre de réception MMSE 


C(f) 


Q'V) 


s q (f)+N 0 n 

La densité spectrale de puissance de la séquence { q (kT h ) } peut se calculer comme suit 


SÆ = jr I \ Qif- f )| 2 

H b k 1 b 

Le récepteur linéaire optimal est donc constitué de deux éléments en cascade 
Un filtre adapté de réponse q (- 1) avec q (t) = g(t) ® h(t) 

Un égaliseur transversal (ligne à délai à cadence 1/ T h dont la réponse est l'inverse de la 
fonction périodique S q (f) + N 0 /2 

La mise en oeuvre de l'égaliseur MMSE optimal implique généralement un nombre infini de 
coefficients (filtre II R). En pratique, on approxime la solution par un nombre fini de coefficients 

{c k }\ = - N C FIR ) 

La conception d’un récepteur MMSE statique ne convient pas en pratique : le canal q(t) n'est pas 
toujours connu (problème de l'estimation de canal), et le canal peut varier dans le temps 

La solution généralement adoptée est la conception d'un récepteur adaptatif qui assure une 
implémentation adaptative à la fois du filtre adapté et de l'égaliseur, de manière combinée. 

'Fractionally-spaced equalizer' : plus facile à mettre en oeuvre, placé directement à la 
sortie du canal, travaille à cadence plus élevée 2 !T h 


Transversal equalizer 


Received 

signal 

<</) 



jKr) 


*— 















Égalisation adaptative et à rétroaction des décisions 

Mise en oeuvre d'un égaliseur adaptatif: requiert le calcul d'une séquence d'erreur e[n ], différence 
entre la sortie courante y[n] de l'égaliseur, et la réponse souhaitée d[n] 


L'algorithme LMS, par exemple, peut être utilisé pour la mise à jour des coefficients w[«] en fonction 
du signal d'erreur e[«] : 


w [ n + 1 ] = w [ n ] + fie [«] x[ri\ 



Output 
‘ vl«] 


Desired 

response 

</(«! 


Jf 


Mode d'entraînement 


au démarrage, une séquence d'entraînement, connue du récepteur, est utilisée par l'émetteur, ce qui 
permet d'initialiser les coefficients de l'égaliseur à des valeurs correctes (disponibilité du signal 
souhaité d\n\ ) 


Mode orienté décisions 

en régime, l'émetteur transmet une séquence inconnue de symboles, et le récepteur prend une 
décision sur les symboles reçus. On peut supposer que ces décisions sont correctes la plupart du 
temps. Les décisions prises par le récepteur constituent alors la réponse souhaitée servant à calculer 
le signal d'erreur. Dans ce mode, le récepteur est capable de suivre les variations lentes du canal. 










































Amélioration de l'égaliseur linéaire MMSE 

Égaliseur à rétroaction des décisions ('decision-feedback equalizer' ou DFE) = égaliseur non-linéaire 
Séquence reçue à la sortie du canal: 

x [/?] = h [0] + X h M a i n - k\ + Z h W a l n ~ k] 

k <0 k< 0 

L'interférence due aux échantillons précurseurs peut être annulée en utilisant les décisions déjà prises 
par le récepteur 



La synthèse de la partie avant ('feedforward') de l'égaliseur a pour but de réduire au minimum (sans 
amplifier le bruit de canal ) l'interférence entre symboles due aux échantillons postcurseurs du canal 

L'interférence due aux échantillons précurseurs sera ensuite parfaitement éliminée par la partie 
arrière ('feedback') de l’égaliseur, à condition que les décisions soient correctes mais l'inconvénient 
de la méthode est le risque de propagation d'erreur 

Un algorithme adaptatif peut être utilisé pour mettre à jour conjointement les coefficients des 2 filtres 
sur base d'un même signal d'erreur 

Pour les canaux fortement dispersifs, l'égaliseur DFE apporte un gain de performance significatif par 
rapport à l'égaliseur linéaire 






















Question 6 : Modulations linéaires 

Décrire et comparer les principaux types de modulations numériques linéaires sans mémoire : PAM, 
QAM, PSK, FSK. Montrer comment ces signaux peuvent être représentés par leurs coordonnées dans un 
espace décrit par des fonctions de base orthonormales, qui constituent un ensemble complet. 
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Un canal est dit sans mémoire lorsque la probabilité d'erreur ne dépend pas de l'instant k d'utilisation 
du canal (k se référant à l’échantillonnage) 

Modulation PAM 

Le principe de la modulation Puise Amplitude Modulation (PAM) réside dans l'encodage du signal 
analogique par rapport à son amplitude. Le signal est échantillonné à intervalles réguliers et le signal 
PAM prendra les valeurs de l'échantillonnage. 

La modulation PAM établit une correspondance entre les symboles Am et les formes d'ondes 
données par 

s m ( t ) = A m g(t) cos (271 f c t) 

pourO<r<r où f c est la fréquence de la porteuse, les symboles notés A m pourm e [1,...,M]. On 
a aussi g(t) une forme d'onde qui influence le spectre occupé par le signal modulé. Le niveau logique 
de la PAM peut varier : 

• PAM-2 : [-1 +1] 

• PAM-4 : [-3 -1 +1 +3] 

• (■••) 

Ainsi, les symboles peuvent prendre les valeurs données par : 

A m (t) = (2m - 1 — M) d 

Où 2 d est la distance entre symboles 

La mise en correspondance des mots de k bits avec les symboles A m peut se faire de différentes 
manières. Il est cependant un arrangement préféré en pratique, appelé encodage de Gray, qui permet 
de minimiser l'erreur de transmission étant donné que le passage d'un symbole à un autre dans la 
trame se fait avec un changement d'un seul bit comme le montre l'illustration qui suit (symboles à 2, 4 
et 8 niveaux) : 
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DSB : Le signal modulé défini par l'équation de la PAM avec porteuse est à double bande latérale 
(DSB- Double Sideband), ce qui n'est pas efficace du point de vue de l’occupation du spectre 

SSB : On peut rendre ce signal à bande latérale unique en lui ajoutant sa transformée de Hilbert : 

s m (t) = A m g(t) cos(2nf c t) - A m g(t) cos{2nf c t ) 

pour 0 <t<T , où g(t) est la transformée de Hilbert de g(t) . Dans ce cas, la bande du signal SSB 
(Single Sideband) n’est plus que la moitié de celle occupée par le signal DSB. 



Baseband Signal 



La modulation PAM peut aussi être employée sans porteuse. On parle alors de PAM en bande 
de base et les formes d'onde sont simplement données par 

Sm (0 = A m g(t) 

Dans le cas particulier de symboles binaires +d et d (PAM-2), on a 

S| (0 = ~s 2 (t) 

Ces deux signaux sont dits opposés et ont même énergie s donnée par 

e = f .5j 2 (t) dt = J s 2 2 ( t ) dt 
o o 

Et un coefficient de corrélation p donné par 



o 


Avec la corrélation p 12 entre deux signaux s { ( t ) 
suit : 



et s 0 ( t ) d'énergie Ej et e 2 qui s'exprime comme 














Modulation QAM 


La modulation Quadrature Amplitude Modulation (QAM) ou modulation d'amplitude en quadrature 
(MAQ) est une modulation où l’on effectue de la PAM sur les deux porteuses en quadrature associées 
à une même fréquence c'est-à-dire à la fois sur le cosinus et le sinus de la même porteuse. Les 
formes d'onde sont données par 

Sm(t) = A mc g(t) cos(2nf c t) - A ms g(t) sin(2nf c f) = V m g(t) ms(2nf c t + 0„,) 

Avec 

V m = y]A m / + A m / et Q m =ATMj% 

Ces différentes équations servent à mettre en évidence que cette modulation peut aussi être vue 
comme une modulation combinée de phase et d'amplitude. En combinant Ml niveaux d'amplitude 
avec M2 valeurs de phase, on a M1M2 symboles dans la constellation. Des exemples de 
constellations sont donnés à la figure ci-dessous (encodage de Gray) 
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Modulation PSK 

La modulation Phase Shift Keying (PSK) ou modulation de phase (MDP) est une modulation pour 
laquelle les formes d'ondes sont données par : 

s m {t) = g(i) cos(2nf c t + 

Où 0,„ = 271 ( ” -1) sont les différentes phases possible de la porteuse qui véhicule l'information. On 
ajoute ainsi un incrément de phase allant de 0 à n au signal en bande de base. Ces signaux 
ont tous de même énergie donnée par 

e = J s m 2 (t)dt =0.5 ]g 2 (t)dt = \ 
o o 

L'ensemble des symboles possibles, que l'on appelle constellation est représenté à la figure qui suit 
pour M = 2,4, 8 ainsi que la mise en correspondance avec les mots binaires par encodage de Gray. 


















Modulation FSK 


Une autre quantité sur laquelle l'information peut être imprimée est la fréquence. Une version très 
simple de la modulation de fréquence numérique, mais suffisante pour en étudier les performances 
dans un canal de type BBGA (Bruit Blanc Gaussien Additif) est celle où les formes d'ondes sont 
données par 

s , n (0 = cos(2nf c t + 2nmA f f) 

Les formes d'onde en bande de base s lm (t) sont données par 

Ce type de modulation porte le nom de FSK : Frequency Shift Keying. L'écart de fréquence entre les 
différentes formes d’onde est donné par A f . Toutes les formes d'onde ont la même énergie qui est 
donnée par 

e = J s m 2 ( t ) dt = y ] cos 2 (2nf c t + 2iunA,t) dt 
o o 

Elles ont des coefficients de corrélation p - kr donnés par (k et r étant des nombres entiers allant de 
1 jusqu'au nombre d'états de symbole possible) 

p kr = 4= f s k ( t ) s r ( t ) dt = j \ cos (2nf c t + 2 nk A f -1 ) cos (2 nf c t + 2nrA,- f) dt = y f cos (2 n(k - r) A f t) dt 
o o o 

P kr = 2Mk-l)A f T sin ( 2K ( k ~ r ) A f T ) 

On obtient un sinus cardinal . Le coefficient de corrélation est bien compris entre +1 et -1 mais 
ici ne peut clairement pas descendre jusqu'à -1. Mais on ne veut pas aller jusque 1 car en modulation 
de fréquence, il faut que les deux formes d’ondes codant 2 états différents ne se ressemblent pas du 
tout ! En effet, noyé dans le bruit, il faut les reconnaître. 

Si on fixe k et r, on a plus que A f T comme degré de liberté avec T la durée du symbole. Si on prend 
A f trop grand, on consomme trop de bande passante, tandis que si on on prend A f trop petit, les 
formes d’ondes se ressemblent trop et ne sont plus distinguable dans le bruit. Le meilleur compromis 
est de prendre A f au premier zéro (racine) du sinus cardinal et donc d’avoir le coefficient de 
corrélation égale à 0. 



Question 7 : Détecteur optimal 

Soit un signal numérique de type PAM-M, pouvant prendre la forme si(t), S 2 (t), sM(t). A partir de la 
représentation construite par le démodulateur; c'est-à-dire les coordonnées rk, comment le détecteur 
peut-il décider au mieux du signal qui se cache derrière ces coordonnées rk ? Notant sm le vecteur des 
coordonnées associées au signal sm(t), le choix se portera sur le signal sm(t) tel que l’on maximise p 
[smlr]. Montrer comment ce critère peut mener au principe du détecteur cohérent. Pourquoi ce récepteur 
porte-t-il le nom de 'cohérent' ? 

Livre de référence 282 - 287 

L'objet de cette partie est d'étudier les structures de réception optimales à utiliser pour un 
signal qui est reçu après passage dans un canal du type BBGA (Bruit Blanc Gaussien Additif). Cela 
signifie qu'à la sortie du canal, on retrouve le signal émis, entaché d'un bruit blanc, Gaussien et additif 
(non corrélé avec le signal émis). Voici le modèle qui nous intéresse 


émetteur 


récepteur 


canal 


s m0) 


► r(t) 


n(t) 


A ne pas confondre : 

Un démodulateur est l'organe qui convertit le signal reçu r(t) en une autre représentation ( par 
exemple, les coordonnées du signal reçu dans un espace dont la dimensionnalité est fixée par la 
modulation avec laquelle on travaill e) et donc quelque chose de plus compact 

Le détecteur est l'organe qui décide, à partir de la représentation produite par le démodulateur, quelle 
est la forme d'onde qui a effectivement été transmise. Le détecteur peut parfois commettre des 
erreurs, et l'objet de cette partie du document est d'évaluer les performances des récepteurs (ici 
en PAM) 

Représentation des signaux 

Un signal peut être entièrement décrit si on connaît ses coordonnées dans un espace décrit par des 
fonctions de base orthonormales, qui constituent un ensemble complet. Les fonctions de base dont 
on a besoin dépendent de la modulation en question. Prenons la modulation PAM-M. On note par N la 
dimension de l'espace auquel il faut recourir 

La modulation PAM requiert une seule fonction de base orthonormale (N = 1 ), donnée par 


fit) = ^ g(t) cos(2nf c t) 


Les différentes formes d'onde possibles, s m (t ), ont des coordonnées notées s m valant 



Dès lors 






Sm(f) Smfif) 



Dans le cas où le signal émis n’est pas corrompu par du bruit, le démodulateur peut projeter le signal 
reçu (c'est-à-dire aussi celui émis puisqu'il n'y a pas de bruit) sur les fonctions de base associées à la 
modulation utilisées et le détecteur peut comparer les coordonnées ainsi obtenues avec celles des 
signaux possibles, et décider. 

Le signal reçu est cependant entaché de bruit, ce qui a pour effet de modifier les coordonnées 
observées par rapport à la situation idéale. En projetant le signal reçu Ht) = s m (t ) +n(t) sur les 
fonctions de base, on trouve alors des coordonnées r k données par 

r k = J dt = 

0 0 0 

avec s mk la m ième coordonnées sur la k ième forme d'onde de base. C'est une variable 
déterministe. Le terme n k est quant à lui une variable aléatoire de moyenne nulle (bruit blanc 

gaussien) provoquée par le bruit. Le terme n(t) étant gaussien, n{k) a une distribution gaussienne. Or 

l'intégrale de n (t) , étant une transformation linéaire, donne une distribution gaussienne. Ces 
coordonnées r k sont donc des variables aléatoires gaussienne de moyenne non nulle, la moyenne 
étant de s mk puisque le bruit n ( t ) est centré et de variance donnée par : 

E [H k ~s mk f ] = E[n k >] =NJ 2 

par définition du bruit blanc ayant une densité spectrale de puissance (PSD) constante. Par ailleurs, 
les coordonnées r k sont aussi décorrélées puisque : 

E[(r k - S mk )(r l - s ml )] = E[n k n,] = 8 a 

avec comme résultat dépendant de k et I : 

si k = 1 , alors 8 kl = 1 donc E[( r k - s mk )(ri - s m i )] = N 0 / 2 

si k ! = 1, alors 8 kt = 0 donc E[( r k - s mk )(r, - s ml )] = 0 

A partir de là, on peut dès lors reconstituer le signal reçu r(t) à partir de ses coordonnées comme suit, 

Ht) = I r k f k (t) + n' (t) = 8 k/ 

k= 1 

Attention, f k (t) permet de représenter tout ce qui peut être transmis par l’émetteur, mais pas 
forcément toute la gamme de bruit ! La collection des r k n'est donc pas suffisante pour reconstruire 
Ht). Il y a aussi le résidu de bruit n(t) à retrouver. Mais n(t) ne contient que du bruit et r k contient 
toute l'information permettant de retrouver les bit transmis et donc n(t) on s'en fiche. En effet, on 
peut montrer que : 

E \r k n' (r)] = 0 

De ce fait, on peut considérer que les r k constituent une représentation équivalente à la 
connaissance de r(t), en tout cas pour ce qui est de l'information dont on a besoin pour prendre la 
décision. Autrement dit encore, la connaissance des coordonnées r k est équivalente à celle de r(t) 
pour ce qui est du problème de détection. 


f. 


s , n dt + J n(t)f k (t) dt = s , + n(k) 


Détecteur optimal 



A partir de la représentation construite par le démodulateur, c'est-à-dire les r k , comment le détecteur 
peut-il décider au mieux du signal qui se cache derrière ces coordonnées r k ? Ce qui intéresse 
l'utilisateur, c'est que le récepteur commette le moins d'erreurs de décision possible. On montrera 
plus loin que le critère d’optimisation du détecteur qui correspond à la probabilité d'erreur minimale, 
est le critère du maximum a posteriori. Cela signifie que l'on décide que s m (t) a été transmis si la 
densité de probabilité de s m (t) conditionnée à r est maximale. Notant sm le vecteur des coordonnées 
associées au signal s m ( t ), le choix se porte sur le signal s m ( t ) tel que l'on maximise 
P [s m I ''] 

qui est la probabilité a posteriori de sm. En utilisant la formule de Bayes, on a encore 

r | -I _ P [ r | s m ] p [,?,„] _ p \r\ p [>„,] 

P \ym I r \ p\j -] m 

I p[r\s m \p[s m ] 

m =1 


p [r | s m ] étant la probabilité de vraisemblance calculable donc connue 

p [ s m ] étant la probabilité à priori. Par exemple en QAM-16 p[s1] = ... = p[s16] = 1 donc connue 
p [ r ] étant toujours constante pour tous les m car cette probabilité ne dépend pas de s m 


Et comme l'algorithme consiste à prendre le maximum a posteriori des p [ s m \ r ] , diviser par 
p\r\ ne change rien , donc on peut simplement maximiser a posteriori le numérateur 
uniquement ! De même, p [ s m ] est identique pour tout m, donc on peut le retirer car on ne cherche 
que le maximum. 

Donc pour résumé, dans le cas particulier où les différents symboles sont équiprobables, on aura 
P [ s m ] = m et dans ce cas, le critère se simplifie. Le choix se porte dans ce cas sur le signal s m (t) 
qui donne lieu à la fonction de vraisemblance la plus élevée. La règle de décision porte le nom de 
décision au sens du maximum de vraisemblance. Ce n'est que dans le cas où les différents 
symboles sont équiprobables a priori que le critère du Maximum A Posteriori (MAP) se ramène au 
Maximum de Vraisemblance (MV). 

Ainsi, dans la formule on a plus que la fonction de vraisemblance p\r\s m ] . En particulier, la 
fonction de vraisemblance A ,(s m ) associée au symbole s m (t) est donc donnée par 


P\r I s m ] = X (s m ) 


1 


N 


<«*•>"” exp( ' - t ?, 




Nn 


Considérons le cas binaire (PAM-2) avec des symboles de probabilités a priori qui peuvent être 
différentes. On décidera que .s 1 , (f) a été transmis si 


pl>iM 
P [ s 2 I r ] 


> 


1 soit encore 


p l>il '•] 
p [ s 2 I r ] 


> 


S 


1 


*2 


Par l'expression de la fonction de vraisemblance on a alors 



gj>iM 

p[s 2 \ r ] 


N 

= exp ( - X 

k= 1 
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Si les symboles sont équiprobables et que donc p [s 2 ] = p [s 4 ] , on décidera que s, ( t ) a été émis si 


N 

I 

k= 1 


N n 


N 

> I 

k= 1 


N n 


autrement dit, si la distance euclidienne entre r et s 2 est plus grande que r et s t alors on choisit que 
a été transmis. 

On voit donc que dans le cas où la perturbation introduite par le canal est un BBGA, et les 
symboles sont équiprobables, le critère du MAP revient à minimiser une distance quadratique 
entre les coordonnées du signal reçu et celles associées aux différents signaux susceptibles 
d'être transmis. On en déduit alors que le récepteur optimal effectue la corrélation entre le signal 
reçu r(r)et toutes les hypothèses et, après correction liée à l'énergie E m du signal émis, 

choisit la plus élevée. 

Rappel : p [r \ s m ] a comme donnée r qui représente un vecteur et s m qui représente la forme 
d'onde numéro m. Par exemple, la probabilité p [r \ s 2 ] donne la probabilité qu'on a transmis la 2ème 
forme d'onde. 

Voici le détecteur optimal au sens du maximum a posteriori. On reçoit un signal r(t) qui est une 
forme d'ondes possibles s m (t). On fait la corrélation de du signal reçu r (t) avec toutes les formes 
d'onde s m (t) possible suivi d'échantillonneurs. Si on a des amplitudes différentes possibles, on doit 
tenir compte des termes de correction E/2. 



























Question 8 : Codage de source 

Définir l'entropie d'une source discrète. A quoi sert le codage de source ? Exposer le premier théorème 
de Shannon, relatif au codage de source. Illustrer ce théorème en comparant plusieurs méthodes de 
codage pour une source discrète dont l'alphabet comporte 4 symboles. 

Slide 1.5-1.11 - Livre de référence 574 - 578 

Entropie d'une source discrète 

Une source discrète sans mémoire peut être vue comme une variable aléatoire discrète émettant des 
symboles indépendant selon un alphabet défini avec un probabilité de p k 

La quantité d'information lié au symbole émis peut être calculé en bit selon la probabilité p k 

dfSk) = log 2 (j- k ) 

Cette quantité d'information ne dépend donc que de p k et peut être aussi considérée comme une 
variable aléatoire de valeur dépendante de p k 

L'entropie de la source H(s) est quant à elle la valeur moyenne de l’information produite par la source 
(considérant tous les symboles possibles de l'alphabet) 

H(S) = E [f(s] = K j!p k log 2 i (JL) 

k = 0 n 

On a donc H(S) = 0 si aucune incertitude n'est liée à l'émission de symbole et H(S ) = log 2 (K) si 
tous les symboles sont équiprobable (p k = k) 

Premier théorème de Shannon 

Soit une source discrète sans mémoire d'entropie H(S ), la longueur moyenne des mots code L 
obtenus par un codeur de source sans distorsion est bornée par: 

L >H(S) 

En pratique aucun codeur n'atteindra ce minimum théorique, il sera donc caractérisé par son 
efficacité de codage r| 


Le but étant d'être le plus efficace possible (éviter la redondance, ne pas avoir de même préfixe qu'un 
autre mot code) 

Codage de source 

Réponse à l'objectif de communication efficace en comprimant au maximum la source 
d'information, on s'interroge sur la façon de représenter efficacement les données produites par une 
source discrète, en connaissant les propriétés statistiques de celle-ci 

Les mots codes produits par le codeur de source sont supposés être sous forme binaire, et le code 
doit être sans distorsion, ce qui signifie que la séquence de symboles originale doit pouvoir être 
reconstruite parfaitement à partir de la séquence binaire codée (pouvant être de longueur différentes) 





Exemple 


Symbole source 

Probabilité d'occurence 

Code 1 

Code 2 

Code 3 

s o 

0.5 

0 

0 

0 

s i 

0.25 

1 

10 

01 

S 2 

0.125 

00 

110 

011 

s 3 

0.125 

11 

111 

01111 


Taille moyenne 

- Code I 1*0.5 + 1*0.25 + 2*0.125 + 2*0.125 = 1.25 

- Code II 1*0.5 + 2*0.25 + 3*0.125 + 3*0.125 = 1.75 

- Code III 1*0.5 + 2*0.25 + 3*0.125 + 4*0.125 = 1.875 

La taille moyenne du code 1 est la plus faible mais le code n'est pas correct, en effet, on ne peut pas 

différencier le ,v 0 (0) avec s 2 (00) 


Entropie de la source 

H(S) = lp k log 2 (j-) = 0.5 * log 2 (l) + 0.25 * log 2 ( 4) + 0.125 * log 2 { 8) + 0.125 * / og ,( 8) 

k = o Pk 

= 0.5 + 0.5 + 0.375 + 0.375 
= 1.75 

L'entropie de la source est donc égale à la longueur moyenne du Code II. Selon le second théorème de 
Shannon, cette valeur moyenne de mot code est l'optimum qu'il est possible d'atteindre. 


Cela confirme aussi la non validité du Code I étant donné que 1.25 < 1.75 (va à l'encontre du 
théorème) 





























Question 9 : Codage de canal 

Définir les notions d'information mutuelle et de capacité de canal. A quoi sert le codage de canal ? 
Exposer le second théorème de Shannon, relatif au codage de canal. Appliquer au canal binaire 
symétrique. Donner l'expression de la capacité d'un canal gaussien. Commenter cette expression. 

Slide 2.16 - 2.31 + Livre de référence 584 - 602 

Information mutuelle 


L'incertitude liée à X se mesure par l'entropie H(X). Étant donné que la sortie Y est une version bruitée 
de X, on cherche à connaître l'incertitude de X en observant Y. Pour l'observation de Y, on peut définit 
l'entropie conditionnelle 


H{X | Y) = X I P (xj, y K ) log 2 ( j ) 

H(X | Y) correspond à la quantité d'incertitude qu'il reste sur X après avoir observé Y 

Mais comme l'entropie H(X) représente l’incertitude liée à X avant d'avoir observé Y, on peut calculer 
l'information mutuelle du canal I(X, Y) représentant la part d'incertitude de liée à X qui a été résolue 
par l'observation de Y 


I(X, Y) = H(X) - H(X | Y) 

On peut aussi calculer l'entropie conjointe 

H{X, Y) = X X P (xj, y K ) bg 2 (^V)) 

K = 0 J = 0 J K 


qui mesure quelle quantité d’information est contenue dans un système de deux variables 
aléatoires 


En ensuite faire le lien avec l'information mutuelle 


I(X, Y) = H(X) + H(Y) - H{X, Y) 


Capacité de canal 

La capacité de canal est définie comme l'information mutuelle maximum entre entrée et sortie du 
canal 

C = max, pix y.I(X, Y) en bit/utilisation du canal 

Avec I(X, Y) = X I P(x:, y k ) log 2 ( ^* |A/) l'information mutuelle 

k = 0j = 0 P yk 







Cas particulier du canal binaire symétrique 

On peut démontrer que l'information mutuelle est maximale pour des symboles d'entrées 
équiprobable p(x 0 ) = p(x x ) = 0.5 

La capacité de canal se calcule alors par C = 1 + p log 2 (p) + (1 - p) Iog 2 (1 - p) 

Pour un canal sans bruit (p = 0) : C = 1 bit / utilisation du canal, correspondant à l'entropie de l'entrée. 
Le canal est inutilisable si p = V 2 : C = 0 



Codage de canal 

Réponse à l'objectif de communication fiable en protégeant le système contre les erreurs dues au 
bruit 

Cela implique l'ajout de redondance contrôlée à l'information utile (opération duale de codage de 
source qui consistait à éliminer la redondance intrinsèque de la source) 

Exemple 

Dans les codes en bloc, on ajoute n - k bits de redondance à des messages de taille k. De cette 
manière, on obtient des mots code de taille n. Le taux de codage est défini comme 

r = L <\ 

n ~ 

L'objectif du codage de canal est de faire en sorte que le message de départ puisse être récupéré à 
partir du mot code reçu en sortie du canal bruité, le tout avec une probabilité d'erreur minime 

Second théorème de Shannon - Channel coding theorem 

Objectif : avoir un système de codage de canal tel que la probabilité qu'un bit de message soit erroné 
en sortie soit arbitrairement faible tout en restant efficace 

Soit une source discrète sans mémoire, d'entropie H(S), produisant des symboles toute les T s 
secondes et un canal discret sans mémoire de capacité C utilisé toute les T c secondes. Alors, si 

H(S) < c 

T s ~ T c 

Il existe un système de codage permettant de récupérer la source avec un taux d'erreur arbitrairement 
faible. Le rapport ■£- est le taux critique 

1 C 



Théorème inverse 

Inversément, si 


H(S) 

T s 


> 


Ç_ 

Te 


Alors, il n'est pas possible de transmettre de l’information de la source sur le canal et le récupérer 
avec un taux d'erreur arbitrairement faible 


Le théorème de codage de canal ne fournit qu'une preuve d'existence de bons code mais ne précise 
absolument pas comment construire ces codes. Il ne précise pas non plus quelle sera la valeur 
pratique du taux d'erreur, il se contente de montrer que le taux d'erreur tend vers 0 quand la longueur 
n du code tend vers l’infini 


Application au canal binaire symétrique 

La source émet des symboles équiprobable (H(S) = 1) toutes les T s secondes. Le canal est utilisé 
chaque T c seconde, le taux de code vaut donc 

r = T c /T s 


Pour être dans les conditions du théorème de Shannon, il faut donc satisfaire la relation r < C 

Pour un canal dont les la probabilité de transition p = 1CT 2 , on a C = 0.9192 . Il suffit donc de choisir 
un taux de code r < 0.9192 pour qu’il soit possible de concevoir un système de codage fournissant 
un taux d'erreur arbitrairement faible 


En comparaison, l'utilisation d'un simple code à répétition implique 
le choix d'un taux de codage r très petit pour atteindre un taux 
d'erreur acceptable 






Capacité d'un canal Gaussien 

Dans le cas d'une variable continue X, de densité de probabilité f x (x ), on peut montrer que l’entropie 
H(X) (définie à la base pour des variables aléatoires discrètes) est toujours infinie. On s'intéressera 
donc à l'entropie différentielle h(x) 


h(X) = J f jx) log 2 [y-J- ] dx 


L'entropie différentielle prend des valeurs finies... par contre, elle peut être négative 

Distribution gaussienne (moyenne u. variance o) 

Entropie différentielle : h(X) = 0.5 log 2 (2 n e a 2 ) 

Dépend uniquement de o 2 , pas de p 

Gaussienne = h(X ) plus élevée que n'importe quelle autre distribution de même variance 

Information mutuelle: l( X ; Y) = h(X) - h( X | Y ), satisfaisant aux mêmes propriétés que l'information 
mutuelle définie pour les variables aléatoires discrètes 


Exemple 

Soit un canal gaussien de bande passante limitée à B et de puissance limitée à P. 

La source X(t), de bande passante B, peut être échantillonnée à cadence 2B pour fournir la 
séquence d'échantillons ( X l , X 2 ,... X K ) pendant une durée T avec K = 2BT 

Le signal transmis sur le canal est perturbé par un bruit AWGN de densité spectrale de 
puissance N 0 /2 dont la bande passante est aussi limitée à B 

Les échantillons du signal reçu s’écrivent Y a = Xk + Nk 

Variance des échantillons de bruit N k : a 2 = N 0 B 

Les échantillons du signal reçu Y k sont statistiquement indépendant 

Contrainte sur la puissance du signal transmis E[X k 2 ] = P 

Capacité du canal gaussien 

C = max fxk(x) {I{X k , Y k ) : E[X k *\= P) 

Avec l'information mutuelle 

K x k> Y k) = KY k ) ~ h(Y k \X k ) = h(Y k ) - h(N k ) 

Le maximum de I(X k , Y k ) sera donc obtenu en choisissant la distribution f x (x) qui maximise h(Y k ) 
tout en respectant la contrainte de puissance P (la distribution gaussienne) 



Le signal reçu est donc gaussien de variance P + o 2 


L'entropie différentielle du signal reçu Y k et du bruit N k sont donc 

h(Y k ) = 0.5 log 2 [2ne ( P + o 2 )] et h(N k ) = 0.5 log 2 (2ne o 2 ) 

La capacité du canal gaussien C vaut 

C = 0.5 log 2 (1 + ^;) [bits / utilisation de canal] 

Si on utilise le canal K fois (transmission de K échantillons X k ) sur une durée de T secondes, la 
capacité par unité de temps devient 

C = y 0.5 log 2 ( 1 + ^ ) = B log 2 ( 1 + ) [bits / secondes] 

Il s'agit là de l'un des plus remarquables résultats de la théorie de l'information, en une simple 
formule, on peut comprendre les effets respectifs des 3 caractéristiques principales du système sur la 
capacité du canal (la bande passante B, la puissance moyenne du signal P, et la densité spectrale de 
puissance du signal reçu NO) 

La capacité du canal dépend de la bande passante B de manière linéaire, alors qu’elle ne dépend du 
rapport signal à bruit ^~b 0 ue de manière logarithmique. Pour une variance de bruit donnée, il est 

donc plus facile d'augmenter la capacité du canal en augmentant la bande passante qu'en 
augmentant la puissance de transmission 

Limite fondamentale d'un système de communication Çchannel coding theorem') : 

Si l'on transmet l'information à un débit binaire inférieur ou égal à C bits/s, il est théoriquement 
possible d'atteindre un taux d'erreur arbitrairement faible. Pour atteindre cette limite, le signal 
transmis devra avoir des propriétés statistiques semblables à celles d'un bruit gaussien. 


Supposons un signal idéal où les données sont transmises à un débit binaire R b égal à la capacité C 
du canal. La puissance transmises vaut donc P = E b C où E b représente l'énergie transmise par bit 

Le système idéal est défini par la relation suivante 

C En E i c/5 _ i 

b = log 2 ( 1 + ) ou encore ^ 

Limite de Shannon (BP infinie) 

lim f* = log (2) = 0.693 

B —> co 7V 0 







Question 10 : Codes correcteurs d'erreurs 

Expliquer le fonctionnement d'un code linéaire en bloc. Quels sont les principaux paramètres qui le 
caractérisent ? Expliquer les notions de matrice d’encodage et de décodage, ainsi que la notion de 
syndrome. Que sont les capacités correctrices et détectrices d'erreurs ? Illustrer à l'aide d'un code 
simple (par exemple Hamming). En quoi diffèrent-ils des codes convolutionnels ? 

Livre de référence 632 - 641 

Définition 

Un code correcteur d'erreurs est une technique de codage basée sur la redondance de l'information, 
destinée à détecter, voire corriger, les erreurs qui peuvent être générées par un canal de 
communication peu fiable. La détection d'une erreur dans un message peut sembler naturelle, sa 
location précise et sa correction sont quant à elles moins évidentes. 

Codes linéaires... 

Dans un code linéaire, si c, etc 2 d ésignent les mots codes respectifs des deux suites 
d'information de k bits d l et d 2 , alors à la suite d'information de k bits d l ® d 1 est associé au mot 
code c l ® c 2 où l’opérateur ® désigne l'addition bit à bit modulo 2. 

... en bloc 

Dans un code en bloc, les n éléments binaires des mots codes sont calculés uniquement avec les k 
bits d'information du bloc courant 

d 2 ... d k \ * [Cj c 2 ... c n ] 

Au besoin, le message initial est découpé en paquets, ou blocs, de taille k. 

Pour un code linéaire en bloc, le mot code c s'obtient à partir du mot d'information d par une 
expression matricielle de la forme : 

c = d. G 

où c est un vecteur ligne de taille 1 x n, d est un vecteur de taille 1 x k et G est une matrice de taille k x 
n appelée matrice génératrice du code 

Le taux d'information r d'un code est le nombre r = k / n 

• r < 1 : ajout d'information 

• r = 1 : pas d'information supplémentaire 

• r > 1 : compression de donnée 

La distance de Hamming entre x et y est le nombre d(x; y) de composantes où x et y 
diffèrent 

Le poids de Hamming w de x est le nombre de composantes non nulles de x 
Exemple 

Si x = 10110 et y = 00101 

la distance de Hamming d(x, y) = 4 
Le poids de Hamming w(x) = 3 et w(y) = 2 


Si on a une liste de mot code, on relève les poids minimum w et distances minimum de 
Hamming d 


Ex : 


mot original 

mot 

codé 

00 

00 

000 

01 

01 

101 

10 

10 

110 

11 

11 

011 


Un code est dit t-correcteur si il permet de corriger toute erreur d’au plus t caractères. La condition 
pour qu’un code soit t-correcteur est la suivante : 

d* >2t + 1 avec t = E 


Forme systématique d'un code 

Un code est entièrement déterminé par sa matrice génératrice G 

c = d. G 

Les vecteurs-ligne de G sont eux-mêmes des mots-codes, et sont supposés linéairement 
indépendants. Cela revient à imposer que si d l 4 d 2 alors c = d x G 4 c 2 = d 2 G . Donc des mots 
d'information différents ne peuvent pas être identiquement codés, en d'autres termes, la matrice G 
est de rang k. 

Les combinaisons linéaires sur les lignes de G ainsi que les permutations sur les colonnes de G 
laissent inchangées la probabilité d'erreur par mot mode : la distance (ou le poids) du code est 
inchangé. 

D'après ces propriétés, la matrice génératrice G peut s'écrire sous la forme (Pivot de Gauss) : 

G = {h\ P kx(n-k)] 

où I k est la matrice identité de taille k et la matrice P est appelée la matrice de parité du 
code. Lorsque G est sous sa forme systématique, les mots code commencent par k bits d'information 
et sont complétés par (n - k) bits de redondance. Donc si d = [d x d 2 ... d k ] alors le mot code devient 

c = dG = I d x d 2 ... d k p l p 2 ... p n _ k ] 


Exemple 

Soit le code linéaire (4,2) binaire : { 0000,1101,0110,1011}. Ce code, ayant quatre mots code, 
peut coder quatre messages qui sont les quatre mots binaires de deux bits : d l = 00 , d 2 = 01, d 3 = 
10, d 4 =11. Faisons ce codage en utilisant une matrice génératrice. Une base pour ce code linéaire 
pourrait être : Cj=1101etc 2 = 0110, ce qui donne : 


C 


110 1 
0 110 









Pour sa forme systématique, on choisit une base comme : c 1 = 1011 et c 2 = 0110 ce qui donne 


G 


10 0 1 
0 110 


Ainsi, par exemple un bloc d = [d x d 2 ] est codé en c = dG = [d x d 2 p x p 2 ]. 

Matrice de contrôle de parité 

On considère un code linéaire C de matrice génératrice G. On appelle matrice de contrôle de parité 
une matrice H de dimension (n - k). n de rang (n - k) (rang plein) qui vérifie : 

H. G t = 0 

Caractérisation d'un mot code 

Si c est un bloc, de mot code d = cG alors on doit avoir 

H d T = H (cG) t = H G t c t = ( HG t ) c T = 0 

Ainsi on a la propriété que le bloc d reçu est un mot code (donc jugé non erroné) si et 
seulement si H d T = 0 

Construction d'une matrice de contrôle de parité 

Soit un code linéaire dont la matrice G est mise sous forme systématique, alors la matrice de 
contrôle de parité H a pour expression 

H = [{P T ) (n - k)xk Un-k\ 

Par exemple, 

0 1 1 0 
10 0 1 


G = 


10 0 1 
0 110 


?< H = 


Décodage par syndrome 

Soit un code linéaire de matrice de contrôle de parité H , c un mot envoyé et y = c + b le mot reçu. 
Alors, on a : 

H y T = H(c + b) T = Hc t + Hb' = Hb T 

Ainsi, le vecteur 5 = H y ' ne dépend que du bruit et non pas du mot envoyé. On l'appelle le 
syndrome associé à y. 

Donc, le syndrome se calcule à partir du mot reçu et de la matrice de contrôle de parité, et ne dépend 
que du bruit et non pas du mot code émis. A partir de la connaissance de H et de y, on peut alors 
espérer déterminer le bruit b car s = H y T = Hb 7 . Une fois le bruit connu, c se déduit simplement par 

c = y ® b . 












Exemples 

Code de parité 

le principe est de rajouter un bit de parité (1 ou 0) de tel sorte que le nombre total de 1 sur le mot code 
formé est paire. Le taux d'information est donc de r = ni . 

Information Mot code 

10111010 101110101 - Nombre impair de 1 

10101010 101010100-Nombre pair de 1 

La matrice génératrice 


Code à répétition 

Il s’agit de la manière la plus "naïve" d'aborder le problème : si la transmission d'un mot 
peut-être entachée d'erreur(s), il n'y a qu'à le transmettre directement plusieurs fois ! 

Information Mot code 

00 00 00 00 

01 01 01 01 

Ici on a donc par exemple 3 couple. Ainsi, s'il survient une erreur dans la transmission, elle n’affectera 
pas les trois couples formant le mot du code. Le décodage peut donc se faire par vote majoritaire sur 
les 3 couples de bits. Pour ce code à répétition (cfr tableau) le taux d'information est seulement 33% 
(33% des bits transmis correspondent à de l'information), tandis que 66% ont pour seul but la 
protection de cette information. 


Code à répétition amélioré 


Information 

00 

01 

10 

11 


Mot code 
00 000 
01 101 
10110 
11011 


Supposons que le mot reçu soit y = 01001. Ce mot ne correspond à aucun mot du code, il comporte 
donc une ou plusieurs erreur. L'idée de la correction est de rechercher parmi tous les mots code celui 
est est le "plus proche". Le taux d'information est de r = 2 . 


Ecart mot reçu/mot du code 
y- 00 000 = 01001 

y-01 101 =00100 

y-10 110 = 11111 
y-11 011 =10010 


Poids de l'erreur 
2 
1 
5 
2 


matrice génératrice 






